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The  rate  coefficients  are  calculated  for  trap 
loss  due  to  excited  state  formation  during 
s  -wave  collisions  of  two  atoms  in  a  light 
field  in  a  cold  atomic  gas  near  conditions 
for  formation  of  a  Bose-Einstein  conden¬ 
sate.  Blue  detuning  from  the  allowed 
atomic  resonance  transition  causes  excita¬ 
tion  of  a  replusive  molecular  potential, 
whereas  red  detuning  causes  excitation 
when  the  light  is  tuned  near  a  bound  vibra¬ 
tional  energy  level  of  an  attractive  molecu¬ 
lar  potential.  In  either  case,  when  the  light 
intensity  is  low  and  the  detuning  is  large 
compared  to  the  natural  linewidth  of  the 
atomic  transition,  the  rate  coefficient  for 
the  collisional  loss  rate  is  proportional  to  a 
molecular  Franck-Condon  factor.  A  simple 
reflection  approximation  formula  is  derived 
which  permits  the  rate  coefficient  to  be 
given  analytically  in  either  case.  The 
Franck-Condon  factor  is  equal  to  |^g(/?c)|2/ 
Dc,  where  ^Fg(Rc)  is  the  ground  state 
scattering  wavefunction  at  the  Condon  point 
Rc,  where  the  quasimolecule  is  in  resonance 


with  the  exciting  light,  and  Dc  is  the  slope 
difference  between  ground  and  excited 
potentials  at  Rq.  The  analytic  reflection 
approximation  formula,  as  well  as  a  simple 
phase-amplitude  formula  for  the  intermedi¬ 
ate  range  wavefunction,  give  excellent 
agreement  with  the  results  of  numerical 
quantum  mechanical  calculations.  The  trap 
loss  rates  due  to  binary  collisions  are  com¬ 
parable  to  or  exceed  those  due  to  atomic 
recoil  heating  for  a  wide  range  of  detunings 
to  the  blue  of  atomic  resonance  and  near 
the  peaks  of  photoassociation  resonances 
for  the  case  of  red  detuning. 
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1.  Introduction 

Atomic  collisions  play  a  crucial  role  in  the  phenom¬ 
ena  of  Bose-Einstein  condensation  (BEC)  [1-3].  Elastic 
momentum-transfer  collisions  control  the  rate  of  evapo¬ 
rative  cooling  leading  to  the  high  phase-space  density 
required  for  the  formation  of  a  condensate.  The  ground 
state  scattering  length,  also  a  property  of  elastic  colli¬ 
sions,  controls  the  nonlinear  coupling  parameter  in  the 
nonlinear  Schrodinger  equation  for  the  condensate 
wavefunction  [4].  Inelastic  collisions  among  condensate 
atoms  also  determine  a  lower  bound  to  the  loss  rate  of 
the  condensate  (other  processes  may  contribute  to  loss 
as  well)  and  thereby  help  determine  the  condensate  life¬ 


time.  Therefore,  it  is  crucial  to  understand  these  ground 
state  collisions.  Calculating  their  rates  requires  that  we 
be  able  to  compute  the  ground  state  scattering  wavefunc¬ 
tion.  Although  these  collisions  occur  in  the  absence  of 
light,  photoassociation  spectroscopy  offers  a  powerful 
probe  of  the  ground  state  wavefunction  [5-8].  Light  also 
has  the  prospect  for  being  an  important  probe  of  con¬ 
densate  properties  and  for  manipulating  a  condensate  or 
even  the  collisions  within  a  condensate.  Therefore,  it  is 
crucially  important  to  understand  how  light  affects  the 
atoms  in  a  cold,  dense  atomic  gas  near  the  BEC  regime. 

This  article  will  be  concerned  with  understanding  the 
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rates  of  binary  collisions  in  a  weak,  off-resonant  light 
field  in  a  cold  atomic  gas.  Since  only  the  case  of  a  weak 
radiation  field  is  treated,  the  effect  of  light  on  the  binary 
collision  can  be  described  perturbatively.  The  case  of  a 
strong  saturating  field  will  be  treated  separately  [9],  The 
detuning  is  assumed  to  be  large  enough  compared  to  the 
natural  atomic  linewidth  yA  and  the  density  N  moderate 
enough  that  binary  events  are  sufficiently  isolated  that 
three-body  collisions  can  be  ignored,  as  well  as  collec¬ 
tive  effects  [10],  Certainly,  considering  binary  events  is 
the  first  step  in  understanding  the  effect  of  light  on 
interacting  atoms.  The  object  is  to  provide  a  simple 
framework  with  which  to  understand  these  binary 
events.  They  would  need  to  be  incorporated  at  a  later 
stage  into  a  more  comprehensive  understanding  of  the 
dynamics  and  stability  of  an  atomic  gas  near  BEC. 
Elsewhere,  we  give  a  brief  account  of  the  possibility  of 
using  photoassociation  spectroscopy  to  probe  many- 
body  interactions  in  a  condensate  [11], 

Another  important  reason  for  having  a  simple  model 
for  understanding  these  binary  collisions  in  a  light  field 
is  to  understand  how  spectroscopic  methods  can  be  used 
as  an  experimental  tool  for  probing  the  ground  state 
scattering  wavefunction.  In  particular,  we  have  proposed 
that  high  resolution  photoassociation  spectroscopy 
would  provide  an  especially  sensitive  probe  of  the 
ground  state  potential  and  scattering  properties  [5],  In 
fact,  this  technique  has  been  used  to  place  bounds  on 
scattering  lengths  for  Li  [8],  Rb  [6],  and  Na  [12,  13] 
ground  state  collisions.  Even  tighter  bounds  should  be 
possible  with  improved  experiments.  In  particular,  we 
will  show  below  how  the  spectroscopy  of  very  cold 
gases  might  be  use  to  measure  properties  of  the  ground 
state  wavefunction. 

The  primary  effect  of  light  on  cold  trapped  atoms  is 
to  cause  loss  processes  by  giving  the  atoms  enough 
energy  to  escape  the  trap.  For  example,  a  single  atom 
scatters  light  of  frequency  a>  detuned  Aa=  oj  -  oja  from 
the  atomic  resonance  frequency  coA  at  the  rate. 


/atomic  — 


(l) 


where  Oa  =  V  lrrl/c  dA  is  the  Rabi  frequency  at  laser 
intensity  I  for  the  atomic  transition  with  transition 
dipole  dA,  and  we  assume  |Aa|  »  flA  and  |Aa|  »  yA;  Aa 
is  defined  so  as  to  be  positive  for  blue  detuning  and 
negative  for  red.  Since  each  scattering  event  gives  a 
recoil  momentum  kick  of  fikw  =  f/oj/c  to  the  atom,  and 
since  the  thermal  momentum  of  the  atom  is  comparable 
to,  or  smaller  than,  f/k,„,  this  atomic  light  scattering 
heats  the  atoms,  causing  their  loss  from  the  trap. 


Binary  events  also  cause  loss  of  atoms  from  a  trap, 
but  with  a  rate  which  is  very  different  for  the  cases  of 
red  and  blue  detuning.  Figure  1  illustrates  the  two  cases 
schematically,  using  the  molecular  potential  energy 
curves  for  the  quasimolecule  comprised  of  the  two  col¬ 
liding  atoms.  Blue  detuning  excites  the  ground  state 
atoms  to  a  repulsive  molecular  state,  causing  the  produc¬ 
tion  of  excited  atoms  with  kinetic  energy  Ee  =  Eg  +  fiAA. 
This  heating  has  been  discussed  in  the  context  of  optical 
shielding  [14,15,9],  The  probability  for  this  free-free 
process  is  a  smooth  function  of  detuning.  On  the  other 
hand,  red  detuning  excites  the  ground  state  atoms  to  a 
bound  eigenstate  of  the  attractive  molecular  state.  This 
only  happens  with  high  probability  if  the  light  is  tuned 
close  to  exact  resonance  with  the  energy  Ev  of  the  bound 
vibrational  state  v.  Since  the  decay  of  the  bound  state 
mostly  leads  to  products  that  are  not  trapped,  this  pro¬ 
cess  gives  rise  to  a  trap-loss  photoassociation  spectrum, 
which  is  now  well  studied  for  the  alkali  dimers  [16-21], 
The  probability  for  this  free -bound  process  is  a  sharply 
peaked  function  of  detuning.  We  have  every  reason  to 
expect  that  a  gas  near  BEC  conditions,  or  even  a  conden¬ 
sate  itself,  will  possess  a  detailed  photoassociation  spec¬ 
trum  for  detuning  on  the  red  side  of  a>A. 

When  the  laser  intensity  I  is  sufficiently  low,  the 
probability  for  both  red  and  blue  detuning  loss  processes 
can  be  described  by  a  Franck-Condon  factor  between 
the  ground  and  the  excited  state  radial  wavefunctions. 
For  both  cases,  we  will  give  here  a  very  simple  analytic 
expression  for  the  Franck-Condon  factor  that  is  in  excel¬ 
lent  agreement  with  fully  quantal  calculations.  We  will 
show  that  the  Franck-Condon  factor  is  proportional  to 
|  WgiRc,  F)|2,  the  square  of  the  ground  state  wavefunction 
at  the  Condon  point  Rc,  where  the  quasimolecule  is  in 
exact  resonance  with  the  light.  At  Rc  the  difference 
between  excited  and  ground  potentials  exactly  matches 
//oj.  Our  expression  for  the  binary  collision  rate  exhibits 
the  proper  quantum  threshold  law  behavior  as  E  — >  0.  It 
predicts  that  the  binary  collisional  loss  rate  per  atom, 
"/binary,  also  varies  as  A A  and  at  typical  BEC  densities 
generally  is  comparable  to  or  exceeds  the  atomic  scatter¬ 
ing  rate,  Eq.  (1),  for  a  wide  range  of  blue  detuning  and 
greatly  exceeds  the  atomic  scatering  rate  when  oj  is 
tuned  near  a  photoassociation  resonance. 

We  will  first  describe  the  model  and  the  basic  proper¬ 
ties  of  the  ground  state  wavefunction.  Then  we  will 
discuss  the  binary  collision  rates  for  the  case  of  blue 
detuning  and  the  reflection  approximation  for  the  free- 
free  Franck-Condon  factors.  Then  we  will  describe  the 
modifications  needed  to  describe  the  case  of  red 
detuning  and  free-bound  Franck-Condon  factors. 
Finally,  we  will  consider  some  of  the  implications  for  a 
cold  atomic  gas  near  BEC. 
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dressed  potentials  cross  at  the  Condon  point  Rc,  the 
point  of  quasimolecular  resonance,  where  in  the  un¬ 
dressed  picture  of  Fig.  1,  the  difference  of  upper  and 
lower  potential  curves  matches  the  photon  energy  hu>\ 


ha a  +  VJRc)  -  VJRC)  =  fico 

(3) 

Ve(flc)  ~  Vg(Rc)  =  fiAk. 

(4) 

In  this  paper  we  will  carry  out  purely  model  illustra¬ 
tive  calculations.  For  this  purpose  it  is  sufficient  to  use 
analytic  Lennard-Jones  potentials  for  the  ground  and 
excited  states. 


V,(R)  =  4s, 


(5) 


Fig.  1.  Schematic  diagram  of  quasimolecular  potential  energy 
curves  for  two  atoms  with  an  allowed  S  —>  P  transition,  showing 
free-bound  and  free-free  absorption  for  respective  red  and  blue  detun¬ 
ing  from  the  atomic  resonance  frequency  (uA- 


2.  Model 

The  binary  collision  rates  for  the  two  cases  described 
in  Fig.  1  will  be  expressed  in  terms  of  Franck-Condon 
factors  between  the  ground  and  excited  state  wavefunc- 
tions.  However,  it  is  very  desirable  to  formulate  the 
problem  using  the  field-dressed  molecular  picture  of  a 
collision  in  a  light  field.  This  procedure  is  well- 
documented  elsewhere  [22-27],  even  for  cold  collisions 
[28-30,5,14],  and  need  not  be  described  in  detail  here. 
A  two-state  model  with  a  single  ground  state  and  a 
single  molecular  excited  state  will  be  assumed.  In  the 
weak-field  limit  it  is  straightforward  to  generalize  to 
multichannel  ground  and  excited  states  for  real  alkali 
systems  with  hyperfine  structure.  In  fact,  we  have  now 
developed  computer  codes  that  let  us  include  both 
ground  and  excited  state  molecular  hyperfine  structure 
for  alkali  species  for  optical  transitions  in  the  weak  field 
limit  [12,  13],  The  ground  state  potential  is  VS(R)  and 
the  excited  state  potential  is  VJR )  -  h  Aa,  where  both  Vg 
and  Ve  -A  0  as  R  — >  The  ground  and  excited  states  are 
coupled  optically  by  the  interaction  matrix  element 
ft  (.HR ),  where  the  molecular  Rabi  frequency  is 


i  =  g  or  e  and  n  =  6  for  the  ground  state  van  der  Waals 
potential  and  n  =  3  for  the  excited  state  resonant  dipole 
potential.  The  long  range  variation  is  CtJRu  and  CdR' 
respectively.  We  use  a  reduced  mass  characteristic  of 
sodium.  It  is  difficult  to  show  on  one  graph  the  poten¬ 
tials  over  the  very  large  range  of  energy  and  distance 
involved  in  cold  collisions.  Figure  2  shows  on  a  logarith¬ 
mic  scale  the  magnitude  \V(R)\/kB  (where  kB  is  the 
Boltzmann  constant)  of  the  long  range  potentials  for 
the  ground  and  excited  states  of  the  Na2  molecule,  for 
which  we  take  Ce  =  1500  atomic  units  (e2al,  where 
e  =  electron  charge  and  aQ  =  Bohr  radius  =  0.  0529  nm) 
and  IC3I  =  10  atomic  units  ( e2al ).  Figure  3  shows  on  a 
linear  scale  the  ground  and  attractive  and  repulsive 
excited  state  potentials,  as  well  as  illustrative  wavefunc- 
tions  for  each.  Since  the  excited  potential  is  orders  of 
magnitude  stronger  than  the  ground  state  at  long  range, 
it  is  an  excellent  approximation  to  neglect  the  ground 
state  contribution  in  calculating  the  Condon  point: 


£3] 

Rc3 


Rc 


C, 

-1a 


1/3 


(6) 


Figure  2  also  indicates  the  detunings  associated  with 
each  R ,  taken  as  a  Condon  point.  Note  that  the  ground 
state  interaction  strength  is  not  negligible  at  intermediate 
R,  but  is  stronger  than  1  p,K  whenever  R  is  less  than 
about  300  a0. 


fl(R)  =  V2t tI/c  d(R)  =  b(R)flA  ,  (2) 

where  d(R)  is  the  molecular  transition  dipole,  and 
0  <  b(R)  <  2/V3  relates  the  molecular  and  atomic 
Rabi  frequencies.  The  particular  value  of  b  depends  on 
the  molecular  states  involved.  The  ground  and  excited 


3.  Ground  State  Wavefunction 

3.1  Long  Range  Form 

We  need  the  ground  state  wavefunction  to  use  in  the 
reflection  formula  for  the  Franck-Condon  factors 
we  derive  below.  When  the  collision  energy  E  is  low 
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Fig.  2.  Magnitude  of  long  range  ground  state  van  der  Waals  and 
excited  state  resonant  dipole  interaction  potentials,  shown  using  log 
scales.  The  interaction  strengths  were  taken  to  be  C&  =  1500  atomic 
units  and  C3  =  10  atomic  units  respectively,  characteristic  of  Na  atom 
interactions.  The  scale  on  the  right  hand  side  of  the  figure  gives  the 
detuning  Z\A  corresponding  to  the  Condon  points  associated  with  the 
resonant  dipole  curve. 


R  (a0) 


Fig.  3.  Long  range  ground  and  attractive  and  repulsive  potentials  on 
a  linear  scale  for  the  same  long  range  parameters  as  for  Fig.  2.  Typical 
free  and  bound  excited  state  wavefunctions  for  a  Na  reduced  mass  are 
also  shown,  for  an  exit  channel  with  5  crrT 1  kinetic  energy  and  for  a 
bound  state  with  a  binding  energy  of  4.87  cm-1.  The  2  p.K  ground  state 
wavefunction  is  nearly  linear  between  60  a0  and  200  (Uv 


enough,  only  collisions  with  zero  relative  angular  mo¬ 
mentum  of  nuclear  motion  contribute  to  scattering  cross 
sections.  The  s-wave  ground  state  wavefunction  takes  on 
the  following  well-known  behavior  as  R  — > 


where  the  vT/jV-ttA2/;,  factor  ensures  energy  normal¬ 
ization,  |(^(£')|1f^(£'))|2  =  S(E  -£").  The  asymptotic 
kinetic  energy  is  E  =  f/2k2J2/j.,  where  /x  =  mj 2  is  the 
reduced  mass,  the  deBroglie  wavelength  is  A,  =  2tt/Ic„, 
and  the  velocity  v„  =  h  A'//x.  It  is  sometimes  more 
convenient  to  write  the  energy  normalization  factor  as 


/  2  n  yg  _  2 

\TTfl2kJ  ~  (hv*,)1'2' 


(8) 


If  E  is  small  enough,  the  asymptotic  phase  rjg  =  mr- 
kgAs,  where  As  is  the  scattering  length  and  n  is  the 
number  of  bound  states  supported  by  the  ground  state 
potential;  the  mr  term  can  be  ignored,  since  it  only 
introduces  an  inessential  phase  factor.  The  scattering 
length  is  a  property  of  the  whole  potential  and  for  real 
atoms  is  very  sensitive  to  small  uncertainties  in  the  short 
range,  chemical  bonding,  range  of  R .  However,  the  scat¬ 
tering  length  provides  a  good  parameter  to  characterize 
the  effect  on  the  long  range  wavefunction  of  the  com¬ 
plete  potential. 

Figure  4a  illustrates  typical  ground  state  wavefunc¬ 
tions  (in  their  real  form  without  the  trivial  complex 
phase  factor  e‘Ve)  with  positive,  zero,  and  negative  scat¬ 
tering  lengths.  These  were  calculated  by  making  slight 
changes  in  the  model  potential  parameters.  The  inner 
region  wavefunction,  say  for  R  <  30  a o,  is  characterized 
by  rapid  oscillations  as  the  local  kinetic  energy 
increases  due  to  the  acceleration  by  the  attractive  ground 
state  potential.  The  intermediate  range  wavefunction, 
say  for  80  a,,  <  R  <  A,/4,  is  nearly  linear  in  R  and 
extrapolates  to  an  intercept  near  R  =  As: 


%(R,  E )  «  e'A 


(2 /yV'2  kJLR  -As) 

Vk. 


(9) 


A  much  more  accurate  representation  of  the  ground 
state  wavefunction  in  the  intermediate  range  is  found  by 
correcting  for  the  variation  in  phase  and  amplitude 
due  to  the  nonzero  ground  state  potential.  The  details 
of  this  derivation,  based  on  an  approximate  treatment  of 
the  rigorous  Milne  equation  for  the  phase-amplitude 
form  of  the  wavefunction,  are  given  in  Appendix  A.  The 
result  is  that  a  more  correct  form  of  the  long  range 
wavefunction  is  found  from  Eqs.  (60),  (69),  and  (72) 

WE>  =  '''■(  Jfc)" 


X  a (R)  sin (k^p(R)),  (10) 

where 


sin(kco/?  +  rj?) 

V£ 


(7) 


a(R)  = 


(11) 
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Fig.  4.  (a)  Comparison  of  characteristic  ground  state  near-threshold 

wavefunctions  at  1.4  |jlK  collision  energy.  The  wavefunctions  were 
calculated  using  a  Na  reduced  mass  and  three  different  hypothetical 
model  potentials  which  yield  scattering  lengths  As  =  -  96  a o,  0,  ao,  and 
+  101  ao.  The  intermediate  range  wavefunction  projects  to  a  node  near 
R  =  As.  Of  course,  the  scattering  length  for  any  real  atomic  species  is 
fixed  by  the  actual  interaction  potential  and  can  not  be  varied  (except 
possibly  by  introducing  external  fields),  (b)  Comparison  of  the  scatter¬ 
ing  wavefunction  for  the  potential  for  which  As  =  +  101  ao  with  the  last 
bound  state  wavefunction  (bound  state  binding  energy  is  2.6  mK), 
when  both  are  given  a  common  WKB  normalization  near  the  distance 
of  the  deepest  part  of  the  ground  state  potential. 


P(R)  =  R-AS- 1  (jr)4  R-  (12) 

Equation  (10)  only  applies  in  the  region  R  »  RB  so  that 
the  correction  term  (RB/R)4  «  1,  where 


and  120  a0  for  a  model  Na  potential  with  As  =  +  84.3  a0. 
Such  a  magnitude  of  ,4,  is  near  the  actual  value  for 
doubly  spin-polarized  Na  [31,32,12],  The  actual  wave- 
function  node  at  87.5  a 0  is  shifted  3.2  cia  further  out  than 
As,  as  accurately  predicted  by  Eqs.  (10)  and  (74).  A 
much  larger  shift  of  1 8  a()  in  node  position  is  predicted 
for  87Rb. 


R(a„) 


Fig.  5.  (a)  Comparison  of  the  exact  calculated  ground  state  wave- 

function  (points)  with  the  predictions  of  the  phase  amplitude  approx¬ 
imation,  Eq.  (10),  using  a  Na  reduced  mass  and  a  potential  for  which 
As  =  84.3  ao.  The  actual  node  position  at  R  =  87.  5  ao  is  shifted  to  larger 
R .  The  node  position  is  independent  of  collision  energy,  as  predicted 
by  Eq.  (12). 

3.2  Short  Range  Form 

Before  leaving  this  discussion  of  the  ground  state 
wavefunction,  it  is  useful  to  comment  on  the  nature  of 
the  short  range  wavefunction  as  E  — >  0.  When  the  colli¬ 
sion  energy  E  decreases  towards  zero,  there  will  always 
be  a  characteristic  quantum  threshold  connection  be¬ 
tween  the  asymptotic  and  short  range  wavefunction.  It  is 
possible  to  define  a  characteristic  distance  Rq  where  the 
WKB  connection  between  the  asymptotic  and  short 
range  wavefunctions  strongly  breaks  down  for  collision 
energies  E  <  Eq  [33,  34].  A  necessary  condition  for 
onset  of  the  threshold  law  behavior  is  E  «  Eq.  The 
values  of  RQ  and  Eq  are  determined  by  requiring  that  the 
maximum  in  the  function  c\X(R,  E)ldR  satisfy  the 
following  condition: 


^=(w^)‘/4  d3)  dA(/?g,  Eq)  1 

V  ’  d R  2' 

is  a  constant  independent  of  energy  with  dimensionality 
of  distance.  It  has  a  value  of  42  ao  for  Na  and  77  ao  for 
Rb.  Figure  5  compares  the  exact  wavefunction  and  the 
result  of  Eq.  (10)  in  the  near  linear  region  between  60  ao 
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where  the  local  deBroglie  wavelength  X(R.E)  is  defined 
by  Eq.  (63)  in  the  Appendix  A.  For  all  E  <  Eq  there  is 
a  range  of  distances  around  RQ  where  dA(R,  E)/dR  >  1/2 
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and  the  WKB  approximation  fails.  Using  Eqs.  (17)  and 
(18)  of  [34], 


Rq  = 


201/4 

21/2(7a)'/6 


Rb  =  1.  043 


(15) 


where  a  =  — 


=  7.0419.  Thus,  we  find  that 


Rb  =  0.96 Rq,  independent  of  potential  parameters  (since 
both  have  identical  scalings  with  Cfl  and  /r  ),  and  the 
defined  parameters  RB  and  Rq  turn  out  to  be  fortuitously 
very  close  in  magnitude.  Since  each  of  these  only  de¬ 
fines  a  qualitative  range  of  distance  associated  with  a 
change  in  the  wavefunction,  they  can  be  used  inter- 
changably.  The  WKB  breakdown  region  is  actually  very 
broad  (An  example  for  a  He  mass  is  in  Fig.  3  of 
Ref.  [33]),  and  becomes  broader  as  energy  is  lowered 
below  Eq  =  19  mK  for  Na  or  Eq  =1.5  mK  for  Rb  [34]. 
Therefore,  significant  departure  from  the  WKB  form 
begins  at  distances  well  to  the  right  of  RQ ,  and  Eqs.  (10) 
and  (11)  show  that  the  amplitude  of  the  intermediate 
range  wavefunction  remains  much  closer  to  its  asymp¬ 
totic  amplitude  1  /VU  than  to  the  local  WKB  ampli¬ 
tude  °c  1  /vk  ( R ,  E).  Thus,  the  amplitude  of  the  ground 
state  wavefunction  does  not  appreciably  exhibit  the  ac¬ 
celeration  due  to  the  ground  state  potential  unless  R  is  on 
the  order  of  or  less  than  Rq  ~  RB.  This  is  one  of  the 
reasons  why  a  semiclassical  analysis  fails  in  the  long 
range  region  as  £  — >  0. 

The  form  of  the  wavefunction  in  Eq.  (60)  is  com¬ 
pletely  rigorous  at  all  distances  and  energies,  providing 
that  the  Milne  equation  described  in  Appendix  A  is 
solved  exactly.  References  [33,  34]  show  that  the  ground 
state  wavefunction  for  the  inner  region  R  «  Rq  is  well 
approximated  by  the  form 


'PUR,  £)  =  (A,kU)1'2  e'A 


(IP)' 

W/T-J 


a^KB(R,  E) 


X  sin(/3^KB(R,  E))  (16) 


=  (Ajkoo) 1/2  'Pf™ +  (R ,  E  =  0),  (17) 


function  near  threshold  and  the  wavefuncton  for  the  last 
bound  state,  both  normalized  at  short  range  to  the  same 
WKB  form.  Scattering  wavefunctions  near  threshold,  as 
well  as  the  last  bound  state  wavefunction,  when  given 
this  common  normalization,  are  almost  indistinguish¬ 
able  for  R  <  Rq.  The  physical  reason  for  this  is  that  the 
local  amplitude  and  phase  are  determined  by  the  ground 
state  potential  in  a  region  where  it  is  very  deep  com¬ 
pared  to  the  initial  collision  energy  or  bound  state  bind¬ 
ing  energy. 

Both  the  inner  range  and  outer  range  approximate 
wavefunctions,  Eqs.  (16)  and  (10),  are  proportional  to 
Vko,  in  accordance  with  the  threshold  law  requirement 
that  the  ground  state  wavefunction  for  R  «  \!k,  every¬ 
where  vanishes  as  VU.  This  ensures  that  the  probabil¬ 
ity  for  inelastic  processes  due  to  the  ground  state  en¬ 
trance  channel  vanish  as  fc*  as  E  -4  0.  This  is  the  basic 
threshold  law  for  inelastic  collision  rates  which  must  be 
satisfied  for  light-induced  inelastic  events  for  either  red 
or  blue  detuning. 


4.  Collisions  for  Blue  Detuning 

The  case  of  cold  collisions  in  a  blue  detuned  light  field 
has  been  extensively  studied  recently,  both  experimen¬ 
tally  and  theoretically  [14,15,9,35-41].  Figure  4a  shows 
the  ground  and  excited  state  potentials  for  the  case  of  a 
very  large  detuning  of  5  cm4,  or  300  GHz,  from  atomic 
resonance.  If  the  light  intensity  is  strong  enough  (this 
will  be  defined  more  precisely  below),  the  ground  state 
atoms  approaching  one  another  encounter  an  effective 
repulsive  interaction  near  the  Condon  point  Rc  and  are 
repelled.  This  gives  rise  to  the  phenomenon  of  optical 
shielding  [15,9,35-41],  by  which  the  light  prevents  the 
atoms  from  approaching  closer  than  Rc ,  thereby  reduc¬ 
ing  the  rate  of  ground  state  processes  which  require  the 
atoms  to  be  closer  together  than  Rc.  The  light  also  mod¬ 
ifies  the  ground  state  elastic  scattering  rate  and  scatter¬ 
ing  length,  and  results  in  the  creation  of  hot  excited  state 
atoms  which  share  the  kinetic  energy  released,  /iA A 
[14,15,9].  Since  we  are  considering  the  weak  field  limit 
here,  the  shielding  effect  is  small  and  we  will  concen¬ 
trate  on  the  latter  energy  release  process.  The  event  rate 
coefficient  for  this  process  is 


where  A,  is  an  energy-independent  constant  having  the 
dimension  of  length  and  a™KB(R,  E)  =  1/Vkg  (R ,  E)  is 
the  WKB  amplitude.  Equation  (16)  shows  that 
the  overall  shape  of  the  inner  range  wavefunction, 
ipwKB  +  (r  «  /?e)(  js  independent  of  energy,  whereas  the 
overall  amplitude  decreases  as  (A,k«)1/2  as  E  — >  0.  This 
is  illustrated  in  Fig.  4b,  which  shows  a  scattering  wave- 


Ke(40  =  PUE ,  40 J}  =  PUE ,  4o),  (18) 

where  the  brackets  imply  an  average  over  the  distribu¬ 
tion  of  ground  state  velocities  v„,  and 

PS(E,  40  =  |5eg(£,  40|2  (19) 
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is  the  probability  for  the  event  by  which  two  ground 
state  atoms  collide  in  a  light  field  and  produce  one 
excited  and  one  ground  state  atom,  both  of  which  have 
enough  kinetic  energy  to  escape  any  weak  trapping 
potential. 

When  the  radiative  coupling  is  small,  the  radiative 
distorted  wave  approximation  [24,26,27],  can  be  used 
for  the  S -matrix  element: 

Seg(E ,  Aa)  =  -27 ri  {V-(E  +  /Id A) |  Veg(fl) I  ¥%E)), 

(20) 

where  Veg(R)  =  fifl(R)  from  Eq.  (2).  The  kinetic  energy 
in  the  asymptotic  ground  and  excited  state  channels  are 


E  and  E  +  /ldA  respectively.  The  asymptotic  excited  state 
wavefunction  is: 


s\n(kc  ,R  +  rje)  . 

(21) 


The  ground  and  excited  state  wavefunctions  are  also 
illustrated  in  Fig.  6a.  The  excited  state  wavefunction 
oscillates  much  more  rapidly  than  the  ground  state  one 
because  of  the  much  greater  kinetic  energy  in  the  exit 
channel. 

By  assuming  that  the  radiative  coupling  V,JR)  is 
either  independent  of  Rox  slowly  varying  with  R  near 
Rc,  then  this  term  can  be  removed  from  the  integrand  in 


R  (a0) 


R  (a0) 


Fig.  6.  (a)  Ground  and  repulsive  excited  state  field-dressed  potentials  for  a  blue  detuning  Aa  of 

5  cm-1  =  300  GHz  and  laser  intensity  7=0.  The  crossing  at  R  =  Rc  becomes  an  avoided  crossing 
when  /  #  0.  Ground  and  excited  state  wavefunctions  are  also  shown,  (b)  Similar  figure  for  the  case 
of  a  red  detuning  Aa  of  -4.87  cm-1,  in  resonance  with  the  v  =  84  level  of  the  attractive  excited 
potential. 
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Eq.  (20),  giving  the  radiative  distorted  wave  result  that 

F 

^  ^  =  hide 

(26) 

PIE,  4a)  =  47t2Vc  40|^(£)>|2 

(22) 

where 

=  4t T2VlFeg{E,  4a), 

(23) 

Dc  = 

^(Ve(R)-Vg(R)) 

(27) 

where  V(:  =  Kg! Rc)  and 

EegOE,  4a)  =  K^(E,  40ln+(£)>|2 

(24) 

1! 

II 

IT 

> 

(28) 

is  a  free-free  Franck-Condon  factor.  Although  direct 
evaluation  of  Feg  by  numerical  quadrature  does  not 
converge,  since  both  wavefunctions  oscillate  with  finite 
amplitude  to  R  =  there  are  several  methods  available 
for  evaluating  Feg  numerically  [27],  as  well  as  approxi¬ 
mate  formulas  based  on  the  local  nature  of  the  integrand 
near  the  Condon  point  R<.  Here  we  use  a  two-state  close 
coupled  scattering  calculation  of  the  collision  in  a  weak 
radiation  field  with  an  A1 -independent  Rabi  frequency, 
and  extract  Seg(e)  from  the  asymptotic  field-dressed 
wavefunction  [24-27],  Then  the  Franck-Condon  factor, 
independent  of  the  laser  intensity  chosen  for  the  calcula¬ 
tion,  is  calculated  from: 


is  the  difference  in  potential  slopes  evaluated  at  Rc  and 
vc  =  v,AR(  )  =  \/2(F  -  VJRcJ/fi  is  the  local  velocity  at 
Rc .  The  second  equation  follows  from  neglecting  Vg  in 
comparison  to  Vc.  The  corresponding  semiclassical 
treatment  in  the  field-dressed  picture  of  Fig.  6a  is  given 
by  the  Landau-Zener  (LZ)  approximation  [14,15,9]  for 
p e: 

Ptz(E,  40  =  2  e“Ac  (1  -  e_Ac)  »  2  Ac.  (29) 

The  latter  approximation  applies  when  the  dimension¬ 
less  LZ  adiabaticity  parameter 


Ccg(£,  40  = 


I  SUE,  40p 

4  t r  Vl 


(25) 


Ac  — 


47t2Vc 

hvcDc 


(30) 


where  Stg(E ,  40  is  calculated  using  a  standard  close 
coupling  scattering  code.  This  method  can  be  general¬ 
ized  to  a  multichannel  version  when  there  are  more  than 
two  channels  [27],  and  thus  could  be  used  for  atoms 
with  hyperfine  structure  [12,  13]. 

Figure  7  shows  Fcg  calculated  for  collision  energy 
E/kB  =  2  /jlK  using  a  model  ground  state  potential  with 
the  reduced  mass  for  Na  collisions  and  a  calculated 
scattering  length  of  +84.3  ciq.  The  detuning  range  corre¬ 
sponding  to  R  =  20  a„  to  200  a„  is  8  THz  to  8  GHz.  The 
Franck-Condon  pattern  shows  an  oscillatory  pattern 
versus  Rc  or  detuning. 

5.  Reflection  Approximation 

We  will  now  turn  our  attention  to  developing  approx¬ 
imations  that  will  give  considerable  physical  insight  into 
the  nature  of  light  scattering  at  very  low  collision  energy. 
Let  us  first  look  at  semiclassical  approximations.  We 
will  assume  a  single  Condon  point  Rc.  Thus,  our  analy¬ 
sis  does  not  apply  to  states  with  significant  contributions 
from  two  Condon  points,  such  as  the  low  vibrational 
levels  of  the  0g  state  of  Na2  [12,  13].  A  semiclassical 
approximation  for  the  free-free  Franck-Condon  factor 
Feg(E)  is  given  by  Eq.  (46)  of  the  reference  [42]  (we 
remove  the  uninteresting  Vc  factor): 


is  small  compared  to  unity.  This  condition  also  ensures 
the  low  intensity  regime  of  linear  variation  of  Pe  with 
light  intenstiy.  Using  Eq.  (25)  to  obtain  F£g(E)  from 
Eq.  (29)  gives  exactly  the  same  expression  as  in 
Eq.  (26).  Note  that  both  Eqs.  (26)  and  (30)  have  left  out 
the  phase  factor  sin2(/3g(/?c)  -  AJRc))-  Actually,  the 
phase  factor  has  been  replaced  by  its  average  value,  1/2. 
Note  that  P\z  given  by  Eq.  (29),  even  with  the  phase 
factor,  does  not  satisfy  the  threshold  law  requirement  to 
vary  as  k,  as  E  0.  In  fact,  Fesgc(E)  from  Eq.  (26)  is  in 
very  poor  agreement  with  the  quantum  mechanical 
feg{E)  in  Fig.  7  and  is  not  shown  in  the  figure. 

The  usual  stationary  phase  derivation  of  Eq.  (26)  can 
be  adapted  to  give  the  correct  threshold  law  as  E  -+  0. 
The  details  are  given  in  Appendix  B.  The  result  is  a  key 
result  of  this  paper,  a  simple  reflection  formula: 

Eeg(E,  4a)  =  yj-  \'F+g(Rc,  Ef.  (31) 

The  Franck-Condon  factor  is  directly  proportional  to  the 
square  of  the  ground  state  wavefunction.  It  only  depends 
on  the  excited  state  through  the  slope  term  Dc .  This 
expression  explicitly  satisfies  the  threshold  law,  since  at 
all  distances  R  «  A,,  |  lF„\2  is  proportional  to  h,  when  E 
is  sufficiently  smaller  than  Eq  [33,  34].  When  R  >  RQ, 
the  quantum  wavefunction  is  well-approximated  by 
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Aa=  3  THz  1  THz  300  GHz  100  GHz 


Fig.  7.  Comparison  of  exact  and  approximate  evaluations  of  the  Franck-Condon  factors 
for  different  Condon  points  Re-  The  correspondance  of  detuning  to  Rc  is  indicated  on  the 
top  horizontal  axis.  The  free-free  Franck-Condon  factors  have  dimensionality  of 
(energy)-2,  where  here  energy  is  expressed  in  cm-1  units.  Since  the  same  magnitude  of  the 
C3  coefficient  was  use  for  the  repulsive  and  attractive  curves  in  our  model  calculations,  the 
exact  free-bound  Franck-Condon  factors  are  seen  to  lie  on  nearly  the  same  curve  as  the 
free-free  ones,  when  the  former  are  divided  by  dEv/dv  to  normalize  them  per  unit  energy. 
This  is  as  predicted  by  Eqs.  (31)  and  (41).  The  reflection  approximation  formula  of  Eq. 
(31),  |^g(Rc)|2/Z>c,  is  in  excellent  agreement  with  the  exact  calculations  over  the  whole 
range  shown.  The  intermediate  range  phase-amplitude  formula,  Eq.  (32),  is  good  agree¬ 
ment  for  R  >  70  do.  (a)  logarithmic  scale  (b)  linear  scale  showing  region  near  outer  nodes. 
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Eqs.  (10) -(12),  and  the  intermediate  range  phase-  the  whole  range  of  temperatures  encountered  for 
amplitude  approximation  to  Fpg  is:  trapped  laser  cooled  atoms,  including  the  E  — •>  0  limit. 


Ecg(E,  Aa)  =  h^D  acsin2(fcxpc). 


(32) 


6.  Collisions  for  Red  Detuning 


where  ac  =  a  (Re)  and  pc  =  p(Rc)-  The  basic  difference 
between  this  formula  and  the  Landau-Zener  one, 
Eq.  (26),  other  than  the  phase  factor,  is  the  presence  of 
the  asymptotic  ground  state  velocity  in  the  denominator 
instead  of  the  local  velocity  it-  If  Rc  is  in  the  range 
Rq  «  Rc  «  A»/2  so  that  the  sin  function  is  linear  in  its 
argument,  then 


F%(E,  Aa)  = 


2 palpi  , 

^¥Dck" 


(33) 


The  equivalent  approximation  for  I\  is  found  from 
Eq.  (23), 


PIE,  Aa)  = 


16? rVl 

hv,  I)( 


Ock^pc- 


(34) 


We  can  carry  out  a  very  similar  treatment  for  the  case 
of  the  free -bound  Franck-Condon  factor  for  red  detun¬ 
ing.  Figure  6b  illustrates  the  dressed  potentials  and 
wavefunctions  for  a  red  detuning  case.  The  basic  differ¬ 
ence  from  blue  detuning  is  that  the  attractive  potential 
supports  discrete  bound  states.  When  the  laser  fre¬ 
quency  co  is  tuned  near  the  position  of  a  particular  bound 
level  with  vibrational  quantum  number  v.  the  level  can 
be  excited  and  decay  to  some  product  p,  which  is  no 
longer  trapped.  The  decay  process  is  primarily  sponta¬ 
neous  emission  leading  to  formation  of  molecular  spe¬ 
cies  or  hot  atoms  [28,  42].  In  the  field  dressed  picture  of 
Fig.  6b,  the  bound  state  is  a  scattering  resonance  embed¬ 
ded  in  the  continuum  very  close  to  threshold.  The  rate 
coefficient  for  red  detuning  is  also  given  by  Eq.  (18),  but 
the  5-matrix  element  is  given  by  a  modified  resonant 
scattering  form  described  by  Napolitano  et  al.  [5]: 


Figure  7  also  shows  the  predictions  of  the  FR  and  FPA 
formulas  compared  to  the  exact  results  for  2  p,K  colli¬ 
sions.  The  full  reflection  formula  is  an  excellent  approx¬ 
imation  over  a  wide  range  of  Condon  points,  and  even 
remains  very  good  into  about  20  a(l  where  the  ground 
state  has  accelerated  the  atoms  to  a  velocity  much  higher 
than  their  initial  one.  The  FPA  approximation  formula  is 
also  excellent  at  intermediate  range,  until  R  begins  to 
approach  Rb~Rq.  Figure  8  compares  the  exact  and  FVA 
results  for  other  collision  energies  up  to  2  mK.  We  see 
that  the  reflection  approximation  is  very  good  for  almost 


\SPg(E,  v ,  T0|2 


yP7s  (E.  V.  Aa) 


E-A , 

h 


+  (yv/2)2) 


(35) 


where  Av  is  the  detuning  relative  to  the  position  E,:  of  the 
bound  state.  The  total  decay  rate  of  the  excited  bound 
state  uis  y„  =  yp  +  ys(E,  v,  Aa)  +  yQ,  where  yp  is  the  rate 
by  which  the  bound  state  resonance  decays  to  the  de¬ 
tected  product,  y ,(E,  v.  Aa)  is  the  stimulated  emission 


Aa  =  10  GHz  3  GHz  1  GHz  500  MHz 


Fig.  8.  Comparison  of  the  exact  free-free  Franck  Condon  factors  and  the  predictions  of 
the  phase-amplitude  formula,  Eq.  (32),  for  higher  collision  energies,  Elk b  =  2  pX,  200  pX, 
and  2000  pX  =  2  |xK.  The  breakdown  of  the  reflection  approximation  is  seen  at  the  larger 
Rc  for  2000  pX. 
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rate  back  to  the  ground  state  continuum,  and  y0  is  the 
decay  rate  due  to  any  other  undetected  processes  (such 
as  molecular  predissociation).  The  critical  element  for 
our  purposes  is  the  ys(E,  v ,  Aa)  factor,  which  describes 
how  the  cold  colliding  atoms  couple  to  the  excited  bound 
state.  For  low  light  intensity,  this  factor  is  given  by  the 
Fermi  golden  rule  expression  [5]: 

y.(E,  v ,  dA)  =  ^  |(^e(u)|yeg(^)|^(£))|2,  (36) 

where  'I'Av)  is  the  unit-normalized  bound  state  wave- 
function  of  the  excited  bound  state.  Using  the  same 
factorization  approximation  as  for  Eq.  (23),  this  can  be 
written  in  terms  of  a  bound-free  Franck-Condon  factor. 


Fes(E,  V,  40  =  f“|  n(Rc,  E)\2,  (41) 

where  'l'+AR<:,  E )  can  be  calculated  exactly  or  taken 
from  Eq.  (10))  above.  The  vibrational  spacing  function 
can  be  evaluated  approximately  from  the  discrete  level 
spacing, 

dEy  _  Ey+  1  E„  |  (A/l\ 

Bv~  2  ’  (  J 

or  more  rigorously  from  the  wavefunction  at  some  short 
range  point  Rn: 


BEy 

dv 


irfr 

2  jJL 


kol^ARo,  v)\2  +  -J- 

Kq 


d  XI'JR»,  v) 


d  R 


(43) 


ys{E,  v,  A A)  = 

O  rrr 

-Jv2c\{Vc(v)\Vt{E))\2  (37) 

0  rrr 

=  Fe,(E,  V ,  Aa).  (38) 


The  assumption  is  that  Vc  is  evaluated  at  a  single 
Condon  point;  if  there  is  more  than  one  Condon  point, 
the  interference  between  the  different  amplitudes  from 
each  Condon  point  would  need  to  be  taken  into  account. 

The  derivation  of  the  reflection  approximation 
formula  for  the  bound-free  Franck-Condon  factor  Feg(£, 
v ,  40  follows  immediately  from  the  same  treatment  in 
Appendix  B  that  we  used  in  deriving  the  free-free 
Franck-Condon  factor.  The  essential  difference  is  that 
the  unit-normalized  bound  state  wavefunction  can  also 
be  written  in  energy-normalized  phase-amplitude  form 
just  like  Eq.  (16),  as  is  commonly  done  in  generalized 
multichannel  quantum  defect  theory  [44],  It  is  only 
necessary  to  introduce  the  vibrational  spacing 


9Ey 

dv 


hVy, 


(39) 


where  vv  is  the  vibrational  frequency,  the  number  of 
complete  vibrational  cycles  per  unit  time.  Then  the  unit 
normalized  bound  state  wavefunction  has  the  form  [44] 


/ flF  \  1/2  /  9 II  \  1/2 

'PAR,  v)  =  (^j  ac(R,  v)  sin(4(/C  v)). 

(40) 

where  a  and  [i  can  be  determined  quantum  mechani¬ 
cally  from  Eqs.  (61)  and  (62),  or  alternatively,  from  the 
WKB  form,  Eq.  (64).  The  derivation  in  Appendix  A, 
adapted  for  a  right  hand  turning  point,  carries  through  to 
give 


where  ko  =  kJRi,).  This  equation  immediately  follows 
from  Eq.  (40)  upon  using  the  WKB  forms  in  Eqs. 
(77)-(78).  The  result  depends  only  weakly  on  the  choice 
of  R0  as  long  as  it  is  in  the  classical  region  of  the 
potential  well.  We  use  Eq.  (43)  in  the  numerical  results 
described  below. 

Although  the  free-free  and  free -bound  Franck- 
Condon  factors  have  different  dimensionality,  they  may 
be  compared  by  evaluating  the  free-bound  Franck- 
Condon  factor  per  unit  energy,  or  Feg(v,  AA)(dEv  ldv)~l. 
Figure  7  also  shows  the  results  of  evaluating  this  func¬ 
tion  for  a  number  of  discrete  bound  levels  using  an 
attractive  potential  that  has  the  same  magnitude  of  the 
C3  long  range  potential  coefficient  as  for  the  repulsive 
state  (see  Fig.  3).  Each  bound  level  defines  a  distinct 
Condon  point  nearly  equal  to  the  outer  classical  turning 
point.  The  free-bound  and  free-free  Franck-Condon 
factors  fall  on  almost  exactly  the  same  curve  in  Fig.  7, 
even  when  Rc  is  as  small  as  20  afl.  This  is  expected  from 
the  fact  that  the  free-free  and  free-bound  transitions  in 
the  model  calculation  have  the  same  ground  state  wave- 
function  and  !)(  at  the  same  Condon  points.  The  results 
in  Fig.  7  show  that  the  reflection  approximation  is  just  as 
good  for  bound  states  as  for  free  states. 

Now  we  can  easily  derive  from  Eqs.  (38)  and  (41)  a 
simple  expression  for  the  stimulated  decay  rate  of  the 
excited  bound  state  to  the  ground  continuum,  ys(E,  Aa): 


ys(E,  Aa)  =  =^^  hvv\  Rc ,  £)|2  (44) 

277  Vc  ,  4  2-2/7  , 

=  T  ~D~c  1VvJnh  sllr(fc»Fc)  (45) 


16t72V5 

hv,Df 


acid  pc 


(46) 


=  PAE)v b. 


(47) 
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The  probability  PC(E )  in  Eq.  (47),  defined  by  the  ex¬ 
pression  in  braces  in  Eq.  (46),  is  identical  in  form  to  the 
probability  in  Eq.  (34)  for  a  free-free  transition.  The 
expression  in  Eq.  (47)  thus  has  a  pleasing  physical  inter¬ 
pretation:  the  total  rate  of  decay  out  of  the  bound  state 
is  the  probability  of  decay  during  one  vibration  cycle 
(corresponding  to  a  “complete”  collision  with  incoming 
and  outgoing  parts)  times  the  frequency  of  vibration 
(number  of  cycles  per  second).  The  overall  probability 
for  the  photoassociation  transition,  given  by  the  reso¬ 
nance  scattering  expression  in  Eq.  (35),  explicitly 
exhibits  the  threshold  law  form  due  to  the  ys  factor  in  the 
numerator. 

The  form  in  Eq.  (44),  proportional  to  the  square  of  the 
ground  state  wavefunction,  will  also  give  the  threshold 
law  form  for  p  -  and  d -waves,  etc.  It  also  explicitly  shows 
why  the  photoassociation  spectra  map  out  the  nodal 
structure  of  the  ground  state  wavefunction  in  the  inten¬ 
sity  pattern  versus  vibrational  quantum  number  [28,  17]: 
as  detuning  Aa  changes  from  level  to  level,  the 
Franck-Condon  factor  tracks  the  ground  state  wavefunc¬ 
tion  at  the  changing  Rc .  This  effect  is  expected  to  be 
much  cleaner  in  very  low  temperature  spectra,  where 
only  ,v  -waves  contribute,  than  in  most  of  the  existing 
higher  temperature  experiments,  where  ,y-,  p  -,  d- waves, 
etc.,  all  may  contribute.  The  p-  and  d -waves  do  not 
generally  have  the  same  nodal  structure  in  the  interme¬ 
diate  range  as  the  s-wave  does. 


This  expression  is  independent  of  collision  temperature, 
as  is  should  be.  Using  Eqs.  (2),  (6)  and  (28)  and  the  fact 
that  we  can  write  [43,  42] 


CW3^A(^)3,  (50) 

the  rate  coefficient  can  be  algebraically  transformed  to 
the  form 


2^lue(AA)  = 


3tt 


(51) 


In  Eq.  (50)  f\  is  a  fraction  of  order  unity  depending  on 
the  particular  molecular  states  involved  (In  Hund’s  case 
(a)  it  is  1/2  for  a  2  state  and  1/4  for  a  IT  state;  in  general 
it  will  be  a  function  of  R  that  can  be  worked  out  for  any 
particular  state).  The  following  definitions  are  used  in 
Eq.  (51): 


bc  =  b(Rc )  (52) 

Na  =  yy  (53) 

Aa 


2  RnW 


1  _  _  _  i 12 

Rc  3  \R 


(54) 


7.  Collision  Rates  Near  EEC 

In  order  to  estimate  the  effect  of  binary  collisions  in 
a  cold  dense  atomic  gas,  it  is  necessary  to  compare  the 
rate  for  the  binary  process  with  the  atomic  scattering 
rate  yatomic  in  Eq.  (1).  The  binary  event  rate  is 

ybinary  =  KJtAJN,  (48) 


Recall  that  b(R)  is  defined  by  Eq.  (2).  The  2bifd3 tt 
factor  in  Eq.  (51)  is  a  dimensionless  molecular  physics 
factor  depending  on  the  specific  transitions  involved, 
and  it  has  a  maximum  value  of 

4_  1 

9tt~  T 


where  Kc  is  given  by  Eq.  (18).  For  the  purpose  of  esti¬ 
mating  the  relative  importance  of  atomic  and  binary 
light  scattering  events  in  causing  loss  of  trapped  atoms 
due  to  heating,  it  is  necessary  to  compare  yatomic  and 
2  Tbinary,  since  two  atoms  are  lost  per  binary  event.  For 
the  large  detuning  and  moderate  density  conditions  we 
consider,  the  influence  of  any  collective  effects  on  light 
scattering  is  not  likely  to  alter  significantly  these  basic 
magnitudes.  We  will  base  our  estimates  on  the  interme¬ 
diate  range  formulas  for  R  >  RB,  Eqs.  (34)  and  (46), 
since  these  apply  to  a  wide  range  of  detuning  from  near 
resonance  to  beyond  100  GHz. 

Using  Eqs.  (34)  in  (18),  the  rate  coefficient  for  trap 
loss  events  for  blue  detuning  in  the  intermediate  range  is 

2K?»%Aa)  =  ahpl  (49) 


The  density  A/,  corresponds  to  one  atom  per  cubic  Aa. 
The  density  near  BEC  conditions  in  a  cold  atomic  gas  is 
typically  larger  than  Nx  .  The  factor  fc  has  a  nominal 
order  of  magnitude  unity  and  reflects  the  phase  of  the 
ground  state  wavefunction.  It  will  have  a  node  near  As  if 
,4,  is  positive  and  in  the  intermediate  range. 

Using  Eq.  (51),  we  find 


2ybinary  g  1  N_ 

Tatomic  7  Nx  ‘ 'C ' 


(55) 


This  fraction  is  nearly  independent  of  the  blue  detuning 
Aa  over  a  wide  range,  since/c  ~  1  (except  near  the  node). 
Since  Nx  =  5  X 1012  cm"3  for  Na  and  2X 1012  cm-3  for  Rb, 
we  see  that  the  order  of  magnitude  of  the  light-induced 
binary  collision  loss  rate  is  comparable  to  or  larger  than 
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the  atomic  light  scattering  rate  when  the  density  is  in  the 
range  of  10°  to  1014  cm-3. 

It  is  also  possible  to  get  a  simple  expression  for  the 
rate  coefficient  for  trap  loss  due  to  a  bound  state  photo¬ 
association  resonance  by  using  Eqs.  (47)  and  (35)  in 
Eq.  (18).  The  thermal  average  can  be  greatly  simplified 
at  very  low  T  when  kBT  is  very  small  compared  to  the 
linewidth  yv  of  the  bound  state  photoassociation  reso¬ 
nance  v.  Since  the  numerator  of  the  resonance  scatter¬ 
ing  5-matrix  element,  Eq.  (35),  satisfies  the  threshold 
law  form,  the  energy-dependence  disappears  in  the 
numerator  of  the  resonance  scatttering  rate  coefficient 
expression.  Then  we  are  left  only  with  the  energy  depen¬ 
dence  in  the  E-Av  term  in  the  denominator  of  Eq.  (35). 
But  since  kBT  is  small  compared  to  either  of  the  terms 
in  the  denominator,  the  E  dependence  is  negligible  and 
can  be  dropped.  We  also  assume  y,  =  yp,  that  is,  each 
decay  results  in  a  loss  event,  since  the  stimulated  decay 
rate  is  assumed  to  be  small  (ys «  y„).  Thus,  the  thermal 
average  rate  coefficient  due  to  resonance  v  becomes,  in 
the  T  — »  0  limit. 


2K?\v)  = 


y 

(A,)2  +\yJ2f 


(56) 


=  2£(e1) 


wyP 

(A)2  +  (yJ2f  ’ 


In  Eq.  (57) 


2 41}(A0  = 


32t73V£  2  2 
~hD^acP'c- 


(57) 


(58) 


is  identical  in  form  to  2kt'"c  from  Eq.  (49).  It  is  the 
single-cycle  loss  rate  coefficient,  due  to  passing  the 
crossing  once  in  each  direction.  The  net  rate  is  deter¬ 
mined  by  the  resonance  expression.  We  can  envision  two 
possible  limits.  The  first  is  the  exact  resonance  case 
when  the  laser  is  tuned  to  the  bound  level  so  A  =  0.  In 
this  case,  the  resonance  enhancement  factor  over  the 
single  pass  result  reduces  to  4vvtv,  where  t„  =  1  lyv  is  the 
decay  time  of  level  v.  But  vvt„  is  just  the  number  Nv  of 
complete  cycles  in  one  lifetime.  So  the  maximum  loss 
rate  coefficient  is  2K'":"(v)  =  8  A' 1  Nv.  The  second  limit 
is  when  the  detuning  is  as  far  as  possible  from  reso¬ 
nance.  Assuming  equal  spacing  of  vibrational  levels,  the 
largest  detuning  is  just  half  of  the  vibrational  spacing; 
thus,  the  maximum  detuning  is  just  A  =  ttvv  »  yv.  This 
gives  1/(7 t2Nv)  «  1  for  the  factor  multiplying  2KtJ>  in 
Eq.  (57).  The  minimum  value  the  rate  coefficient  can 
take  is  2 K^m  =  2K[]  ’/(tt2Nv).  Of  course,  we  must  add  the 
off-resonance  contributions  of  all  levels  v—l,  v  +  1, 
etc.,  to  get  the  total  of /-resonant  rate  coefficient.  If  the 
vibrational  spacing  is  assumed  to  be  uniform,  the 


minimum  rate  coefficient  from  summimg  over  all  such 
contributions  is  =  2K[1,/(41V„).  Thus,  we  estimate 


4^-  (2^)  <  (2£/d)  <  (4A1,)(2^C1)).  (59) 


From  Eq.  (59)  and  (55),  and  the  fact  that  Nv  »  1 
whenever  Aa  »  yA,  we  see  that  2y^ary  will  be  a  sharply 
peaked  function  that  is  much  smaller  than  2yb]“ary  (for 
the  same  |z\A|)  when  at  is  tuned  between  resonances,  but 
which  is  much  larger  than  2 7^! Ay  when  the  co  is  tuned 
to  a  photoassociation  resonance.  Therefore,  by  detuning 
to  the  red  side  of  atomic  resonance,  the  effect  of  binary 
collisions  can  be  made  to  be  either  very  large  or  very 
small  compared  to  yai,.„nc-  The  cold  atomic  gas  will  have 
a  prominent  photoassociation  spectrum,  and  the  rate  for 
photodestruction  at  the  peak  positions  will  measure  the 
ground  state  pair  correlation  function  °c|'/'gCRc)|2  at  the 
Condon  points  for  the  transitions.  There  is  a  discrete  set 
of  Condon  points  //  Ju  ),  corresponding  to  the  outer  turn¬ 
ing  points  for  the  vibrational  levels  of  the  upper  molec¬ 
ular  state.  There  are  several  well-known  alkali  dimer 
transitions  that  have  already  been  studied  [7]  in  much 
hotter  traps  than  the  <  1  p,K  traps  that  are  now  available. 
The  intriguing  question  remains  whether  3-body  or  col¬ 
lective  effects  might  perturb  the  binary  photoassociation 
spectrum,  causing  broadening  or  shifts  in  the  resonance 
line  profiles.  Kagan  et  al.  [45]  have  suggested  the  possi¬ 
bility  of  observing  quantum  correlations  in  optical  prop¬ 
erties  of  low  temperature  spin-polarized  gases.  The 
photoassociation  spectrum  of  a  cold  dense  atomic  gas 
near  BEC  might  be  a  very  fruitful  source  of  experimen¬ 
tal  and  theoretical  investigation  of  2-  and  3-body  interac¬ 
tions  in  such  systems  [11]. 

Certainly  the  analysis  in  this  paper  shows  that  light 
can  induce  significant  loss  processes,  in  comparison  to 
that  caused  by  free  atom  light  scattering,  due  to  colli¬ 
sions  between  atoms  in  a  cold  atomic  gas  near  1  p,K  or 
below  at  densities  >  1013  cm-3.  If  the  density  should 
increase  to  >  1015  cm-3,  the  collisional  loss  rate  will 
dominate  that  due  to  atomic  recoil  heating,  except 
between  photoassociation  resonances  for  large  enough 
red  detuning.  For  the  case  of  the  MIT  Na  trap  [3]  the 
blue  detuning  of  the  laser  used  to  plug  the  hole  in  the 
magnetic  trap  was  so  far  off-resonance  that  any  atomic 
or  binary  light  scattering  did  not  prevent  the  formation 
of  a  condensate.  But  for  much  smaller  detunings  the 
simple  Franck-Condon  formulas  we  have  derived  here 
give  a  way  to  estimate  the  binary  collision  rates  for  trap 
loss  via  excited  state  production.  More  detailed  treat¬ 
ment,  including  incorporating  the  actual  quasimolecular 
structure  of  alkali  dimer  states,  could  be  readily  incor¬ 
porated  into  this  framework.  We  have  also  developed  a 
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theory  for  treating  the  nonperturbative  effect  of  light, 
including  multiphoton  effects  and  the  modification  of 
the  ground  state  scattering  length  [9],  Careful  attention 
to  collective  effects  also  needs  to  be  given  for  the  case 
of  higher  densities.  In  any  case,  it  seems  very  unlikely 
that  a  condensate  could  be  stable  for  very  long  in  the 
presence  of  light  tuned  within  a  few  hundred  GHz  of 
atomic  resonance. 


Therefore,  we  want  to  calculate  the  departure  of  this 
large,  nearly  independent  of  R ,  amplitude  from  its 
asymptotic  limit.  We  do  this  by  approximating  the  last 
term  of  Eq.  (61)  by  1/a3  »  ki  a.  Then  Eq.  (61)  becomes 

^+(k2(R)-kJ)  a  =  0  (66) 


8.  Appendix  A.  Milne  Equation  Analysis 
of  Asymptotic  Wavefunction 

The  use  of  the  exact  Milne  quantum  mechanical  equa¬ 
tions  for  the  phase  and  amplitude  of  the  wavefunction 
lets  us  do  an  exact  analysis  of  the  asymptotic  wavefunc¬ 
tion.  Let  the  wavefunction  be  writen  in  phase-amplitude 
energy-normalized  form: 


d2a  2ll  Tr  ,ns  . 

dR~-~~^  V(R)a  =  0  (67) 

d 2  a  ( 2p,  \  1 

~AR2  =  ~  \Jf  CV  R*a'  (68) 

Equation  (68)  can  be  immediately  integrated  to  give 


V(R,  E)  = 


(2p_Y2 

wr-J 


a(R ,  E)  sin(/3 (R,  £)), 


(60) 


1  - 


2 pC6  \ 

20  fi2R4) 


(69) 


where  a  satisfies  the  Milne  equation  [46,  47], 

+  (6I) 
The  phase  f3  is  calculated  from  a  by 

«fi-£)  =  /s(SLE7<ifi-  (<i2) 

and  the  local  wave  vector  is 

k(R,E)=  yJ^(E-V(R)).  (63) 

The  asymptotic  potential  is  V(R)  =  -  G/Rfi. 

In  the  semiclassical  picture,  we  have  the  identifica¬ 
tion  of  a  with  the  WKB  amplitude 


(70) 


where  RB  is  defined  by  Eq.  (13).  When  R  =  201,4Rb  = 
2.1  Rb,  the  amplitude  a ( R )  has  only  changed  by  5  % 
from  its  asymptotic  value  a*.  When  R  =  RB.  the  correc¬ 
tion  term  is  unity,  and  the  asymptotic  approximation 
used  in  setting  up  the  approximate  Eq.  (66)  breaks 
down.  Therefore,  the  amplitude  expression  in  Eq.  (70) 
can  only  be  applied  for  R  »  RB. 

Now  we  can  turn  our  attention  to  deriving  the  correc¬ 
tion  to  the  asymptotic  phase.  This  is  done  by  substituting 
the  above  amplitude,  Eq.  (70),  into  the  phase  [Eq.  (62)]. 
The  asymptotic  phase  at  some  suitable  large  value  of 
R  =  R  ,  (a  distance  well  beyond  the  range  of  the  poten¬ 
tial)  is  /T  =  kx(R*,-As),  where  As  is  the  s-wave  scatter- 
inglength.  The  phase  function  at  R  is  calculated  by  inte¬ 
grating  backwards  from  R„, 


aWKB(R ,  E )  = 


1 

Vk(R,  E )  ’ 


(64) 


P(R)  =  kAR,-A,)  + 


a(Rf 


d  R 


(71) 


and  the  above  forms  are  equivalent  to  the  WKB  wave- 
function.  The  exact  solution  to  the  quantum  equations 
[Eqs.  (61)  and  (62)]  give  an  exact  quantum  wavefunc¬ 
tion.  Here  we  will  use  an  asymptotic  analysis  for  the 
wavefunction  at  moderately  long  range. 

At  long  range  the  amplitude  becomes  constant,  and 
very  large  as  £  0: 


=  k„(Ro°-As) 


d  R 


~  koo(R„-As) 


1  +  2 


d  R 


a„ 


1 

Vg 


(65) 


=  UR-As-l  ^)R 


3  \R 


(72) 
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Combining  the  corrected  phase  and  amplitude 
expressions,  Eqs.  (72)  and  (70),  in  Eq.  (60),  we  obtain 
the  expression  given  in  Eq.  (10)  in  the  text.  Note  that  the 
normalization  stays  close  to  its  asymptotic  value,  in  spite 
of  the  fact  that  the  local  k(R.E)  can  be  much  larger  than 
h,  in  that  part  of  the  moderately  long  range  region 
where  |Vg|  »  E  (see  Fig.  2).  Thus,  the  semiclassical 
local  WKB  normalization  does  not  apply.  Note  also  that 
the  node  of  the  wavefunction  at  R„  near  A,  for  the  posi¬ 
tive  scattering  length  case  is  shifted  from  its  value 
projected  from  the  asymptotic  As.  This  is  because  the 
phase  near  As  has  not  yet  reached  its  asymptotic  limiting 
value.  This  shift  in  the  actual  node  to  larger  R  can 
readily  be  estimated  from  Eq.  (72): 

Rn-As-\(f\R„  =  0  (73) 

and  R„  ~  As  implies  Rn  ~  A.,^1  +  ^  .  (74) 


Feg(E,  Aa)  =  \leg(E,  Aa)|2,  (79) 

where 

7eg(£,  A A)  =  I  agaesin/3gsin/3ed R  (80) 

Jo 

=  ^  J  agac  cos ((3g  -  /3e)d7C  (81) 

In  the  last  line  we  have  made  the  standard  approximation 
of  neglecting  the  rapidly  oscillating  phase  sum  term. 
The  phase  difference  is  expanded  in  a  Taylor  series 
about  a  point  Rp: 

<KR,  E  )  =  p g(R,  E)  -  p C(R,  Ec)  (82) 

=  b0  +  b,(R  -RP)  +  \  b2(R  -  Rp )2  +  .  .  .  ,  (83) 

where 


9.  Appendix  B.  Reflection  Formula 

This  derivation  of  the  reflection  approximation  ex¬ 
pression  follows  Jablonski  [48]  and  the  Appendix  of  Ref. 
[49],  We  begin  by  writing  the  ground  and  excited  state 
wavefunctions  in  energy-normalized  phase-amplitude 
form 


%(R,  E)  = 


(  2al\  1/2 


ag(R,  E)  sin(/3g(R,  E)) 


(75) 


b0  =  Pg(RP.  E)  -  fiARp,  Ee) 

(84) 

=  &0 RP,  E)--^ 

(85) 

d/3 g(R.  E)  d Pe(R.  Ee) 

(86) 

d  R  rp  d  R  Rp 

II 

R 

"5  t— ‘ 

5 

1 

?T- 

bi 

(87) 

~  k„-  ke(Rp,  Ee) 

(88) 

'PAR.  Ec)  = 


(  1/2 

\7T^2/ 


aJR.  Ee)  sin(/3e(R,  Ee )). 


(76) 


b2  = 


d2/3 AR.  E ) 
d  R2 


d 2 PAR,  Ee) 


d/C 


(89) 


where  Ee  =  E  +  //Aa.  The  ground  state  amplitude  and 
phase  are  given  by  Eqs.  (70)  and  (72),  whereas  the 
excited  state  quantities  can  be  taken  to  be  the  standard 
WKB  ones: 


d ke(R,  Ee) 

d  R  Rp 

(90) 

pA), 

PkeiRp)  ' 

(91) 

aJR.  Ee)  = 


1 

VkAR.  Ee) 


(77) 


PAR.  Ee)  =  f  UR',  EA  d R'  +  ^7r,  (78) 

where  R,  is  the  inner  turning  point  of  the  repulsive 
potential.  The  Franck-Condon  factor  is 


In  the  normal  semiclassical  theory,  Rp  is  taken  to  be  the 
Condon  point,  where  kg  =  kc  so  that  b:  =  0.  Here  we 
simply  take  the  point  where  b i  vanishes.  Since  the 
excited  state  potential  is  so  steep  with  respect  to  the 
effectively  flat  ground  state  potential,  it  is  an  excellent 
approximation  to  take  Rp  =  R,  =  Rc  to  all  be  equivalent. 
Then  from  Eq.  (78),  /3e(Rp,  E)  =  tt/4.  The  expression  for 
b2  follows  from  neglecting  the  ground  state  contribution 
in  comparison  to  the  excited  state  one. 
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We  thus  evaluate 

/eg (E)  =  ag(Rc,  E)ae(R ,  Ee)  |  J  cos (b0  +  y  x2)dx 

/  Trfi  ^  \  1/2 

=  (  YjjTh)  E)  sin(/3g(/?c,  £)).  (92) 

where  we  have  used  J  cos^b0  ±  -y-  dr 

=  ^[COS(Z?0±?)-  (93) 

Using  Eqs.  (75)  and  (92)  in  Eq.  (79)  gives  the  desired 
reflection  approximation  to  the  Franck-Condon  factor: 

Feg(E,AA)=^-\%(Rc,E)\2.  (94) 

Of 
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